Dirty Hubbard rings: renormalization group and exact diagonalization studies 
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We studied mesoscopic Hubbard rings with impurities using the renormalization group (RG) 
technique and the exact diagonalizations. The exact diagonalization calculations showed that the 
charge stiffness normalized by the value of clean system has a peak as a function of electron-electron 
interactions. Previous RG analysis has succeeded to show the enhancement of the stiffness at weak 
interactions, but does not show the peak behavior. We derived RG equations with 4&f impurity 
scattering term, which was ignored in the previous studies, and reproduced the overall behavior of 
the exact diagonalization results. 



I. INTRODUCTION 



Low dimensional systems often show interesting quantum phenomena, and there particle-particle interactions and 
jy^ [ impurities sometimes play crucial roles. Mesoscopic rings are known to be one of the good examples and have been 
studied over the last few years in particular. When an external flux is applied through a mesoscopic ring, the current 
flows along the ring and the magnitude and the direction of the current changes as a function of the applied flux. The 
current flows in equilibrium and is therefore called persistent current. The presence of persistent current was first 
I ' predicted by Byers and Yang, and Biittiker et al. Q, and observed in several recent experiments |Q. The magnitude 
, of the observed persistent current is, however, much larger than the theoretical value calculated for non-interacting 
electrons with impurities |3), and therefore we need to take account of other factors. The factors, other than impurities, 
that might affect the persistent current are, for example, particle-particle interaction, periodic potential, spin degree 
of freedom, and multichannels. If the system has only one of these, the situation is rather simple and we believe we 
already well understand it. Impurities would reduce the persistent current ||]. Interactions would not change the 
persistent current [^,^). Periodic potentials would change the shape of the persistent current as a function of the 
applied flux, but would not change its magnitude very much ||. Spin degree of freedom or multichannels just give 
some numerical factor to the magnitude of the persistent current. 

Now let us consider the system that has two or more of the above factors. The role of each factor could be different in 
this case. In the presence of particle-particle interactions and impurities, for example, the persistent current decreases 
as the interactions increase , while the interaction is an irrelevant parameter when it exists with no other factors. 
In general [EJ, interaction becomes relevant if it coexists with translationally noninvariant factors. As this example 
shows, the situations when two or more factors coexist are usually complicated and therefore interesting. In this paper 
we focus on the system with impurities, interactions, periodic potentials, and spins. The most typical and maybe 
the simplest model that contains all of them would be dirty Hubbard ring. Dirty Hubbard ring has been recently 
studied analytically |?],[| and numerically (HJjyJ . Giamarchi and Shastry H showed using the renormalization group 
(RG) analysis that repulsive interactions between particles would enhance the persistent current while they would 
j_j ■ suppress the persistent current in spinless fermion systems fl2|-|l4|l. The persistent current does not continue to be 
enhanced, however, as we increase the interactions further. The exact diagonalization calculations Jl^,|ll| show that 
the persistent current has a peak at a characteristic interaction strength and begins to be suppressed by the interaction 
beyond that point. Two origins are expected to give the suppression. One is periodic potential. Even in the absence 
of impurities, interactions would suppress the persistent current on a lattice and the suppression is more emphasized 
as the particle density approaches to the half filling. This effect would exist in the dirty systems. The second origin 
that would suppress the persistent current would be the Akp Fourier component of impurity scatterings. For strong 
electron-electron interactions, Hubbard model approaches to a spinless fermion model, and the Akp scatterings in 
Hubbard model corresponds to the 2kp scatterings in the spinless fermion model. Since it is known [[12-14 that the 



2Rf scatterings in the spinless fermion model would suppress the persistent current, we expect the Akp scatterings 
would also suppress the persistent current of Hubbard model. The previous RG equations do not contain the AkF 
scattering term and are not complete to reproduce the numerical results. We will derive RG equations taking account 
of the AkF scatterings and will show the peak behavior of the persistent current. 

The paper is organized as follows. In Sec. 2 we show the results of our exact diagonalization calculations. RG 
analysis is given in Sec. 3. Section 4 is devoted to summary of this study. 
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II. NUMERICAL RESULTS 



We first show our exact diagonalization results, which are basically the same calculations as the previous ones 
pll| , pl| . We used the modified Lanzcos method |ll| to calculate the ground state energy E of Hubbard ring: 

H = -t e l2 ^ 4,oL clc ja + c.c. + Uj2^n lL +J2v imp {i)n l<7 (1) 

<ij>,a i i<7 

as a function of the applied flux (f>. We assumed the uniform probability for the impurity potential Vi mp between 
— W and W . The impurity average of the energy E was taken over 512 to 3072 realizations. The persistent current 
I(4>) is given by the first derivative of the averaged energy with respect to the flux cf>, and the charge stiffness D c is 
proprtional to the second derivative of the energy. Since it it easier to calculate the charge stiffness than to calculate 
the persistent current in the RG analysis, and since the charge stiffness is a measure of the persistent current j7|,pT[, 
we here show the results of the charge stiffness calculations. In the calculations we always chose even number of 
particles per spin. The reason is the following. For simplicity, let us consider free electrons with the numbers of up 
and down spins, N-\ and N±, being equal. The ground state energy of the free electrons changes as we change the 
particle number per spin, N a , from even to odd, and it has the relation E even (<f>) = E dd(4> + <^o/2) where 4>o is the 
flux quanta. If we define the charge stiffness by D c = (L/2)d 2 E((f>)/d(f) 2 \ c f >= o, the stiffness as a function of the particle 
density, does not change smoothly when N a moves between even and odd. We can avoid this problem if we take the 
second derivative at 4> = when N a is odd, and at <j> = <po/2 when N a is even. This causes another problem, however. 
In U — > oo limit Hubbard model becomes the free spinless fermion model, where the particle number is always even 
since iVj = iVj,. Therefore we should take the derivative of the ground state energy at <f> = 4>q/2 in U — > oo limit 
although the derivative has to be taken at <fi — at small U. If we choose N a to be even, the system is always in the 
'even number' sector in both limits of U — and U — > oo, and we do not have to change the definition of the stiffness. 
For this reason, we chose even number of particles per spin and used the definition D c = (L/2)d 2 E((f>)/d(j) 2 \ < p =r f >0 /2 in 
the following exact diagonalization calculations. 

Figure ([!]) shows the charge stiffness D c as a function of U/t in the system of L=6 and N-\ = N± = 2. As observed in 
the previous exact diagonalization studies Jl(],[ll| , it has a peak, and the peak position shifts towards higher U jt for the 
stronger impurity potential. The enhancement at small U is due to the impurity scatterings, since we know the charge 
stiffness in the clean systems, D C (W — 0), decreases monotonically as U becomes larger, If we normalize the stiffness 
by D C (W = 0), the peak still exists and the peak positions are located almost in the same place for the different 
impurity strengths (Fig. (Q)). Since the normalization effectively eliminates the periodic potential contribution, Fig. 
(||) suggests that the reduction of the charge stiffness at large U is not only due to the periodic potentials but also 
due to the impurities. In the next section we will try to reproduce this behavior of the charge stiffness using the RG 
technique. 



III. RG CALCULATIONS 



One dimensional Hubbard model with impurities can be written in the boson representation by the following 
Lagrangian: 



where 



L = L p + L„ + L g + L 2 k F + L^kp, (2) 



L v = ^ K J dx{u y 1 {d T 4> v {x,T)) 2 +u u (d x ^ u (x,T)) 2 }, 

On r 

L 9 = /r, ^2 / dxC08(V8<j><r(x,T)), 

{2-Kay J 



L 



2 k 



— f dx^ 2 (x)e iV2 'l'^ x ' TS > +2kFX cos(V2<£ Jx,r))+h.c, 
na J 

L ikF = — [ dx^(x)e iV5 ^ (x ^ +ikFX + h.c, 
na J 

v = p,a, and a is the lattice constant. The last two terms in Eq. (|J) represent the 2kp and 4kp Fourier components 
of impurity scatterings, respectively. The parameter g\± is equal to U a, where U is the on-site repulsion of Hubbard 
model. We assume gaussian distribution for the random fields £2 and £4: P(^) = ^2 exp(— J \^2{x)\ 2 dx /W2) and 
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P[^i) = A r 4exp(— J \^4{x)\ 2 dx/W4), where N 2 and N4 are the normalization factors. We did not include the forward 
scattering due to impurities because it can be eliminated by redefining the field 4> p ||. 

With help of the replica trick, we can integrate out the random fields M and the final form of the action is 



S — S p + S a + Sg + S 2 k F + S4k F , (3) 



where 



^=E^/ dxdr{u-\d T 4>iix,r)f+u v {dMx,T)f}, 

i 

5 s=E(t^I / ^cos(V8^(x,t)), 



S 2fcF ^r^ / d^Tdr'cos(\/2^(x,T))cos(\/2^(x,rO)cos(\/2(0;(x,T)-^(x,r')), 

ij 

Sik F = ^2 ( Wi \2 I dxdrdr' cos(V8((j) p (x , r) - <^,(x, t )), 
ij 

where i and j are the replica indices. Now we are ready to derive the RG equations. Up to the lowest order of gi±, W 2 
and W4, the RG equations are obtained in the following form: 



dK p (l) _ 1 K z p u p 



dl 2 \ u a 

du p (l) _ 1 / K p u 2 p 
dl 2 u„ 



(A 2 (0 + A 4 (0), 
(A 2 (0 + A 4 (0), 



^M = - l -KM 2 (l)+y{in 

-^>=--K a u a A 2 (l), (4) 

dy{l) 

dl 
A 2 (Q 

dl 
A 4 (0 



dl 



{2-2K a {l))y{l)-W 2 {l), 
(3-K p (l)-K a (l)-y(l))A 2 (l), 
(3-4if p (0)A 4 (0, 



where y = gi±/jTu cr , A 2 = (2W 2 a / '■nv? a )(u a j 'u p ) K '» , and A4 = %W^au a j-Kv? p . If we ignore the 4fcp scatterings, i.e. if 
we set W4 = 0, the RG equations reduce to the ones derived by Giamarch and Schulz 

We solved the Bethe ansatz equations of clean Hubbard ring [|l6|,|l7j to determine the initial values for the integration 
of the RG equations. From the Bethe ansatz equations we obtain the stiffness and the susceptibility of charge and 
spin sectors. The parameters K's and u's can then be calculated from the following relations 

K P = ^^2DpY P , 



2 Dp_ 



K a = W2AxX CT , (5) 

4 [d7 

where D v and \v are i respectively, the stiffness and the susceptibility of the v(— p, a) sector and p is the particle 
number density. Note that the first two equations in Eq. (^) are exact, while the last two are approximate and valid 
when U/t is small 0,| . 
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Using the obtained parameters as the initial value of the RG equations, we integrated the equations numerically 
and stopped when I reached log(L/a). Figure (||) shows the normalized charge stiffness D C /D C (W — 0) calculated 
with W4 = 0. As obtained by Giamarchi and Shastry ffi], the stiffness grows as the interactions increase, and never 
drops down. When W4 is finite, however, the normalized charge stiffness has a peak as seen in Fig. The peak 
stays almost in the same position for the different impurity potential strengths, and this feature is also observed 
in the exact diagonalization calculations. We therefore believe the RG equations (^) are complete to reproduce the 
qualitative feature of the exact diagonalization calculations. 

IV. SUMMARY 

The previous RG calculations that contain the 2kp impurity scatterings succeeded to show the enhancement of 
the charge stiffness by the interaction U, which is an opposite feature to the spinless fermion systems where the 
interactions would suppress the stiffness. As we increase U further, however, the exact diagonalization studies show 
the stiffness starts to decrease making a peak at a characteristic U. The peak behavior still exists even when we plot 
the normalized stiffness D C (W)/ D C (W = 0) as a function of U. This feature cannot be obtained from the previous 
RG equations and we need to take account of the missing term in their calculations. 

We derived the RG equations for the dirty Hubbard rings with the Akp as well as the 2kp impurity scatterings being 
included, and calculated the charge stiffness D c for various impurity strengths W and electron-electron interactions 
U. At U = 00 the Hubbard model is equivalent to the free spinless fermion model, where fermi wave vector is twice 
as large as the one of the original Hubbard model. Therefore the 4fc f scatterings become dominant as we approach 
U/t — » 00 limit. Since the RG and the numerical calculations for the spinless fermions showed interactions would 
suppress the charge stiffness in dirty systems, we can expect that the RG equations with the 4fcp term should give the 
decreasing behavior of the charge stiffness with increasing U /t. We numerically integrated the RG equations and the 
normalized stiffness of the results makes a peak at a certain U/t, while it grows monotonically without the 4&f term. 
It is also found that the peak position does not change very much for different impurity strengths. These results are 
consistent to the numerical analysis. 

If we ignore the 4fc^ term in the RG analysis, the persistent current becomes close to the impurity- free case at large 
U . In real systems, however, the 4Jtp scatterings do exist and the persistent current cannot reach the impurity-free 
value. Therefore, one-dimensional Hubbard rings may not be sufficient to explain the large persistent current observed 
in the experiments , and we should take account of other factors such as multichannels. 

Let us close this paper by commenting on the Coulomb interaction cases. Since arbitrary weak Coulomb interaction 
in one-dimensional electron systems makes the 4fc^ density correlation dominant and drives the system to a Wigner 
crystal |fis| , the 4fci? component of impurity scatterings always plays central roles, and the persistent current and the 
charge stiffness would always be suppressed by the Coulomb interactions as a consequence. The exact diagonalization 
calculations |lf[] actually show that the Coulomb interactions would suppress the persistent current unless the impurity 
potentials are too strong |lS| ]. 
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FIG. 1. Charge stiffness vs U/t for W = 0.1,0.5, and 1.0 from top to bottom. The system size is L — 8. The particle 
numbers are Nf = N± —2. 



FIG. 2. Normalized charge stiffness vs U/t. The parameters are the same as in Fig. (1) 



FIG. 3. Normalized charge stiffness vs U/t when A4 = 0. The system size and the particle number per spin are L = 100 
and Nf = Ni = 22 respectively. The initial value of A2 is set to be 1, 30, 60 (xlO -4 ) from top to bottom. 

FIG. 4. The normalized charge stiffness vs U/t when A4 = 4A2. The system size, the particle number, and A2 are the 
same as in Fig. (3) 
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